USITP 2001-8 
solv-int/0110021 



A unified treatment of quartic invariants at fixed and arbitrary 

energy. 

Max Karlovini* 
Department of Physics - Stockholm University 
Giuseppe Pucacco^ 
Dipartimento di Fisica - Universita di Roma "Tor Vergata" 
Kjell Rosquist^ and Lars Samuelsson^ 
Department of Physics - Stockholm University 



Abstract 

Two-dimensional Hamiltonian systems admitting second invariants which are quartic in the momenta 
are investigated using the Jacobi geometrization of the dynamics. This approach allows for a unified 
treatment of invariants at both arbitrary and fixed energy. In the differential geometric picture, the 
quartic invariant corresponds to the existence of a fourth rank Killing tensor. Expressing the Jacobi 
metric in terms of a Kahler potential, the integrability condition for the existence of the Killing tensor 
at fixed energy is a non-linear equation involving the Kahler potential. At arbitrary energy, further 
conditions must be imposed which lead to an overdetermined system with isolated solutions. We obtain 
several new integrable and superintegrable systems in addition to all previously known examples. 



PACS numbers: 03.20.-fi 11.30.-j 



*e-mail: max@physto.se 
te-mail: pucacco(S)roma2. infn.it 
ie-mail: kr@physto.se 
§e-mail: larsam@physto.se 



1 



1 Introduction 



The study of Hamiltonian systems by means of geometric techniques provides fruitful clues about their 
integrability ^. In particular, geometrization of the dynamics by using the Jacobi metric [Q, ||], is 

a standard tool to turn a natural Hamiltonian system into a geodesic flow over a suitable Riemannian 
manifold. Therefore, it is natural to investigate integrability of two-dimensional systems by looking for 
invariants corresponding to Killing tensors of a conformal Riemannian geometry. These geometric objects in 
fact directly produce invariants polynomial in the momenta: the degree of this polynomial is the same as the 
rank of the Killing tensor. Since the energy enters in the Jacobi metric as a parameter, the conditions for 
the existence of a Killing tensor may happen to be satisfied only at a fixed value of energy. In this case we 
speak of integrability at fixed energy in distinction to the usual notion of integrability which is understood 
to be valid at arbitrary energy. 

In [P, quadratic invariants at arbitrary and fixed energy (respectively, strongly and weakly conserved 
phase-space functions) for 2-dimensional Hamiltonian systems were treated in a unified way. The integrability 
condition for quadratic invariants involves an arbitrary analytic function S{z). For invariants at arbitrary 
energy, the function S{z) is a second degree polynomial with real second derivative and the integrability 
condition then reduces to the classical Darboux's condition for quadratic invariants at arbitrary energy [^. 
Thereafter, the possibility of searching for linear and quadratic invariants at fixed energy was addressed and 
some examples of system admitting a second quadratic invariant at zero energy were provided. 

Weak invariants, also called configurational invariants, have been discussed by Hall and by Sarlet, 
Leach and Cantrijn js). Hietarinta ||^, in his account of the direct methods for the search of the second 
invariant, also provides a review of all the known systems admitting one or more configurational invariants. 
Generalizing the approach of paper , Karlovini and Rosquist JlO[ have discussed the existence of invariants 
cubic in the momenta at both fixed and arbitrary energy. Beside giving a list of all known systems admitting 
a cubic strong invariant, they also find a superintegrable system admitting a cubic invariant related to an 
energy dependent linear invariant. In the present paper we are going to discuss the case of the quartic 
invariant associated with the existence of a fourth-rank Killing tensor. Fourth-rank Killing tensors have 
previously been studied by one of the authors (see [|ll|) in the Lorentzian case. 

In analogy to the third-rank case, it turns out that it is natural to introduce a Kahler potential for the 
Jacobi metric. Expressing the conditions imposed on the geometry by Killing tensor equations in terms 
of this Kahler potential, there remains a system of nonlinear partial differential equations. Imposing the 
condition of strong integrability, it turns out that in general this system becomes overdetermined, so that, like 
for cubic invariants, only isolated mechanical systems exist with a second invariant quartic in the momenta. 
Due to the complexity of the system of equations, it is impossible to get a fully general solution. However, 
with our approach we are able to find some new integrable and superintegrable cases. In particular, one may 
wonder if the knowledge of a weak invariant can provide information about the global dynamical behaviour 
of the system. We introduce a method which, starting from a family of weakly integrable systems, leads to 
finding a higher-order invariant for isolated members of the family which are therefore strongly integrable. 
The second invariant depends on the energy as a parameter. 

The plan of the paper is as follows: in section 2 we recall the necessary and sufficient conditions for 
the existence of a Killing tensor in the Jacobi geometry and their application in problems of analytical 
mechanics; in section 3 we examine the particular case of fourth-rank Killing tensors; in section 4 we analize 
the conditions for strong integrability with a quartic invariant; in section 5 we present the method by which, 
starting from a family of weakly integrable systems, isolated examples of strongly integrable systems with 
quartic second invariant can be found; in section 6 we give tables of integrable systems admitting a quartic 
invariant which include new integrable and superintegrable systems in addition to all known cases; in section 
7 we give our conclusions. 

2 Jacobi geometry and Killing tensors 
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2.1 Geometric representation of the dynamics 

We are interested in the classical two-dimensional systems with Hamiltonian function 

n^\{pl+pl) + V{x,y). (1) 

The approach followed in studying the integrability properties of these systems is based on the Jacobi 
geometrization procedure, which associates to the Hamiltonian flow produced by (|l|), a geodesic flow on a 
Riemannian manifold endowed with a positive definite metric 

ds''j=g^0dq"dq'', a,P=l,2. (2) 

It can be shown (see, e.g. |^), that the Jacobi metric 

gaf3^2(E-V)hc,f3, (3) 

where hap is the metric of the flat space, generates a geodesic flow corresponding to the natural mechanical 
system with potential V = V{q) at energy E. In fact, introducing the geodesic, or "Jacobi", Hamiltonian 

(4) 



the geodesic equations 

dg" dUj 1 



Pf3, 



dsj dp^ 2{E-V) 
dpa _ dUj _ 1 dV ^' 

'dTJ ~ ^Ikt ~ ~2{E-V) dq"' 
are equivalent to the canonical equations of motion given by the Hamiltonian (|^), 

dq" _ on 
dt dpa ' 

dpa _ _&H ^ ' 

IT ~ ~dq^' 

if the natural identification of canonical coordinates 

q^^X, q^^y, Pl^Px, P2^Py (7) 

is made. The phase-space trajectories are parametrized by the affine parameter sj that is related to the 
standard time variable by 

dsj ^2[E -V{q{t))]dt. (8) 
The Jacobi geometry corresponding to Hamiltonian (|l]) is then 

ds^ = 2G{dx^ + dy^) = 2Gdzdz, G = E-V, (9) 

where we have introduced null variables z — x + iy and z — x ~ iy. 

All tensor calculations will be done in the standard null frame defined as 



= ^/Gdz, QP ^ VCdz, (10) 

in which the metric takes the simplest possible form 

ds^ = 2dn°dn°. (11) 

We use the convention that tensor indices in this frame take the values and 0, while in any coordinate 
frame the values will be the names of the coordinates (e.g. z and z). 
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2.2 Invariants polynomial in the momenta and Killing tensors 



One of the standard tools of the geometric approach is the investigation of integrabihty by looking for 
invariants generated by Killing tensors. Let us see the simple case of second-order Killing tensor equations 

^(a/3;7) = 0. (12) 

The existence of a second-rank Killing tensor, that is a symmetric tensor which satisfies eq.(p^, implies the 
existence of a conserved quadratic function 

/,/ = K^^PcPfS (13) 

commuting with the Jacobi Hamiltonian. In fact, it is easy to check that the vanishing of the Poisson 
bracket of this function with the Hamiltonian implies eq. ([T^) . In analogy with the theory of second order 
invariants, the Jacobi geometry approach leads, as it stands, to the determination of higher-order invariants. 
Generalizing eq.(Q2), an m'th rank Killing tensor is a symmetric tensor satisfying equation 



-f^(AxiM2...p™;M™+i) = 0- (14) 
It gives rise to the m 'th order (in the momenta) invariant 

Ij ^ K''it^--^^"^p^,...p^^. (15) 

Two remarks are in order here: 

1) The function ( p^ is a weak invariant, in the sense that, in general, it is a conserved quantity in the 
dynamics fixed by the given value of the energy appearing in the definition of the Jacobi geometry (|^) . To 
grant it the full title of strong invariant, it must satisfy the additional requirement of satisfying the Killing 
tensor equations at arbitrary energy. 

2) The function (p^), let us call it the Jacobi invariant, is a homogeneous polynomial of degree equal 
to the rank of the corresponding Killing tensor. To transform it into the ordinary invariant in the physical 
time gauge there is a straightforward recipe consisting in replacing the parameter E appearing in // with 
the corresponding Hamiltonian function (|l]). As a consequence, the physical invariant 

^ = ^j\E^n (16) 

becomes a polynomial which is either even or odd in the momenta. 

Karlovini and Rosquist (2000, |jl^) have discussed the existence of third-rank Killing tensors, giving a 
list of all known integrable or superintegrable systems admitting a cubic strong invariant. Here we are going 
to discuss the case of the quartic invariant associated with the existence of a fourth-rank Killing tensor. 
Moreover, we want to exploit the results concerning the existence of higher-order invariants to generate new 
solutions starting from the assumption of the existence of a second-rank Killing tensor. 



3 Quartic invariants corresponding to fourth rank Killing tensors 

In this section we derive the necessary and sufficient intcgrability condition for the Jacobi metric to admit 
a fourth-rank Killing tensor at a fixed value of the energy parameter E. In the following section we proceed 
by finding the conditions that ensure that the Killing tensor equations be satisfied at arbitrary values of the 
energy, with the restriction to the case in which the energy dependence of the Killing tensor is such that the 
corresponding invariant of the physical Hamiltonian is quartic in the momenta. As usual in our approach 
(see ||ll[), we decompose the fourth-rank Killing tensor in the form: 

Kaf:!jS = PafljS + P{af39y5) + gi^afSg-fS) (17) 

or equivalently 

KafjfS = PcxPiS + Kj^ap9-,S) " Z^g(af}gi8) (18) 
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where Kap — K'^ap-y and K — K'^-y. Pap-yS is the trace-free (conformal) part, Pap — P'^ap-y its trace and 
K the fuU trace of the KiUing tensor. In this way, the KiUing tensor equations are spht into the trace-free 
"components" and the equation for the trace 

P{aP-fS:fi) — T^P" (ap-i-u9Sti) = 0; (19) 



P(ap;-i) - -j^P" (a-ug-yS) + -j^P" aPr^ = 0, (20) 

K^a = -^P^c-p- (21) 

It is advantageous to employ also the coordinate frame components of the conformal part and of its trace 
when parametrizing the five independent components of KapjS- Thus, the components of the Killing tensor 
can be written as _ 

^0000 — Lr — Lr i> 

^0000 = ^^00 = Igp'' = \gW^ 

Koom = \k (22) 
^0000 = 2"^™ ^ -GP^-^ = -GR{z,z) 

T:fzzzz cf 

^0000 — — Lt O, 



where, as usual, S = S{z) is a holomorphic function, so that S = S{z) comes from (19), but the function 
R = R{z, z) is a generic function of the arguments. With this parametrization, the equation for the trace 
( ^ ) has the integrability condition 

2(G,,,i? - GrzzR) + m,zR.z - G^Az) + G{R,,, - R,,) = 0, (23) 

which is analogous to that obtained for the second rank case (but now R is not holomorphic!). 
We now proceed to solve the Killing tensor equations. Eq.(EOh gives 



S-^/^R,, + 2(5i/^G),, = (24) 

and its complex conjugate. In analogy with the approach followed in |l^, we make a coordinate transfor- 
mation to put the conformal Killing tensor in the simplest (standard) form. Using a conformal transformation 
of the form 

w = H{z), 

w = X + iY, (25) 

z = X + iy, 

the transformation of the conformal tensor is such that 

S{w) := = [H'{z)]^P''''' = [H'{z)]*S{z). (26) 

If we make the standard choice 

S{w) := P"""""" = 1, (27) 
the conformal transformation is then provided by the function 

H'{z) = [S{z)]-^/\ (28) 

Eq.(^) becomes 

(5-i/2i?),,a + 2 =0, (29) 
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with 

G^\H'{z)\^^G^ ^/\S\G. (30) 
Moreover, under the conformal transformation generated by eg. (p8|) , the function R transforms as 

R = P"""" = [H'{z)]^P'' S-^/^R. (31) 

Inserting 7? into eq.(p9|), we have the set 

R.w + 2 G w — 0, 

(32) 

R,w + 2 G,u, — 0. 

The solution of this system can be found in terms of a real Kahler potential IC(w,iv): 

G = /C,^«. (33) 

R 2 JCy^WW- 

whereas the integrability condition (p3| ) determines the following equation for the Kahler potential: 

^ ,www^ ,ww ^" .www^ .ww^w^ — 0- (34) 

This is the necessary and sufficient condition for the existence of a fourth rank Killing tensor. In analogy 
to the third rank case (see [0) and in contrast to the first and second rank case Q, the condition is highly 
nonlinear. The same condition has been already found in and treated in ||], even if not in the context of 
the present geometric approach. Note that in eq.(7.5.15), there is a misprint for a factor 2 missing. 
In the standardized coordinate frame the second invariant can be written in the form 

I J = 2 ^{pi + 2Rnjpl} + ^n^K, (35) 
whereas in the original null coordinate frame the second invariant is 

Ij = 2^{S pt + 2 RHj pl} + ^njK, (36) 
where ^ ^ ^ 

is the Jacobi Hamiltonian expressed in the three different reference frames. The trace function K is found 
by integrating the system 

^i^s, + i?,„G + 2i?G,^ = 0, 

o _ _ (38) 

-K^^ + R^G + 2RG^u, = 0- 

4 Arbitrary energy invariants 

Eq. ( |3^ ) gives the the necessary and sufficient condition for a two-dimensional Riemannian geometry to admit 
a fourth-rank Killing tensor. If we interpret the geometry according to eq.(^) as that providing a natural 
mechanical system, we can derive the additional conditions that make the Killing tensor equations satisfied 
for every value of the energy E. We use the following form for the Kahler potential, 

/C = £;[zz + 2 3f?{A(z)}]-1', (39) 

where A is a holomorphic function independent of E and the real pre-potential is such that 

= V. (40) 
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In fact, no loss of generality is implied by letting the energy dependence of the Kahler potential be 
prescribed by eq.(|39|) and by taking the function S4{z) to be energy independent, given that we only take 
interest in the cases for which the physical invariant / is a quartic polynomial in the momenta. Since the 
analogous statement was made without proving it for the cubic case we give a proof here. We begin by 
noting that the physical invariant /, obtained from the Jacobi invariant Ij of eq. ( ^ ) according to prescription 
(p^), can be written as 

I^2?fi{Spt + Bplp,}+C{p,p,f, (41) 

where 

B = 2RG-\ 

3 (42) 

Since the energy parameter E is now assumed to have been replaced by the physical Hamiltonian Ti, the 
functions S, B and G will be dependent on the momenta and pz, but clearly only through the combination 
PzPz that appears in 7i. This implies that all five terms on the right hand side of eq.(p]) have to be quartic 
polynomials in the momenta since / is quartic by definition. This can be clearly seen by viewing / as a 
function of the two independent momenta functions PzPz and pz/pzi rather than the momenta Pz and Pz 
themselves. Indeed, / can be recast into the form 

2 

i={pzP,f yQk{pz/p,)\ (43) 



k=-2 



where 



Q-2 — Q2 ^ S, 

Q-i=Q~i = B, (44) 
Qq = C. 

Clearly, with the coefficients of [pz/Pz)^ only depending on the momenta through PzPz, there can be no 
cancellation of possible non-quartic polynomial dependence in the individual terms. Hence, we conclude 
that the functions Spz, Bp'^pz and G{pzPz)'^, as well as the complex conjugates of the first two, must be 
quartic momenta polynomials. It immediately follows that S, B and C are polynomials in {pzPz)~^j with 
coefficients depending only on z and z, of degree zero, one and two, respectively. This directly proves the 
part of statement about S having no energy dependence. Moreover, since G — E — V turns into T — 2pzPz 
when E is replaced by H = T + V, it now follows from eqs.(^2|) that the functions R and K are restricted 
to be first respectively second degree polynomials in PzPz, which is the same as saying that they must be 
the same type of polynomials in the energy parameter E, before making the substitution E Ti. Using 
the second of eqs.(|3^) and the conformal transformation formulae given by eqs.(p8|) and (^ij), the function 
R can be expressed in terms of S and /C according to 

R = -2{SIC,zz + \s'JC.z), (45) 

while eqs.(^8[), which notably are form invariant under conformal transformations, give that the trace K is 
given by integrating 

^K^z + RzG + 2RG^z = (46) 
and its complex conjugate. Now, there is clearly no loss in generality to write the Kahler potential /C as 

K,^ E[zz + 2^{K{E,z)}]~'^{z,z), (47) 

with A being analytic in z and 5* being energy independent. Moreover, since S is energy independent and 
i? is a first degree polynomial in E, it follows from eq.(^5|) that we are able to write A(i?, z) as 

K{E,z)=k^{z)+A^{E,z), (48) 
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where Ai(i?, z) does not contribute to eq.(^5|). However, since eq.(^6|) shows that A(z) enters into K g only 
through R, we might as well set Ai{E,z) to zero, which proves the fact that we could have taken A{E,z) 
to be energy independent from the outset. We may finally note that with both R and G being first degree 
polynomials in E, we find from eq.(^6|) that Kg is a second degree polynomial in E, which indeed is a 
necessary and sufficient condition for if to be a second degree polynomial as well, up to addition of an 
irrelevant integration constant which in principle can have an arbitrary energy dependence. This completes 
the proof that the form of the Kahler potential given by eq. ( |39| ) is the most general one which is needed for 
full generality. 

We are now in the position to solve the integrability condition at arbitrary energy. In the original z 
coordinates, eq.(p^) takes the form 



Substituting the ansatz (39) into eq.(|49|), this condition provides a second degree polynomial in E: 

A2E^ + AiE + Ao = 0. (50) 

The coefficients Aq, Ai, A2 must vanish separately if the equation must hold for arbitrary energy. Therefore 
we obtain the three equations: 

A2 = q[ [A'"S + ^A"S' + \iz + A')S"U} = 0, (51) 

+ ^S'{{z + A')*,...- + ^(*,.. + A"*,,,-)) (52) 
+ i5"(*,, + (z- + A')*..)] J = 0, 

Ao - 3{ + 2^',,,,-^',,,) + i5'(*,,,j^',, + ^^..j^',..) + ^S"^^,g^,,] = 0. (53) 

It turns useful to express the system of eqs.(^-^3|) also in the transformed coordinates since this simplifies 
the computations when the S function is of higher degree. Using the notation 

F{w) = H-\z{w)) (54) 

for the inverse conformal transformation, the coefficients Aq, Ai, A2 become 

A2 = 3{ [F"'FF'F' + 2{F"fFF' + A'" F' F' + A" F" F'] (55) 

Ai = 3{ [{F"'F + A"')*,^fi, + 2{F"F + A")^ + F'F'^^^^^ + 2F"i?'*,^„] J^, (56) 
Applying the differential operator 



dzdz 

to eq. (pl|) , we get the condition 

3{5""(z)} = 0, (58) 
so that the form allowed to the S function to have integrability at arbitrary energy is 

S{z)^ az^ + I3z^ +-iz'^ + 5z + e, a e R, /3,7,(5,eeC (59) 
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In 1^ it has been shown that, to get arbitrary energy quadratic invariants in systems with two degrees of 
freedom, the function S{z) must satisfy the condition 

= 0, (60) 

that is it must be a second degree polynomial with real second derivative. Together with the result ( |5^ ) and 
the corresponding one obtained in the third-rank case ||lO|| , 

^{S'^\z)} = 0, (61) 

we can guess that, as a general rule, the analytic function S{z) representing the conformal part of a Killing 
tensor of arbitrary rank m is required to satisfy the condition 

5R|(-^)"5(z)| =0 (modd), 

dz (g2) 

)™5(z)| =0 (meven). 

As a simple illustration of the solution of the above set of eqs.(^Tl)-(|5^), let us take 

A = 0, S ^b, beR. 
Eq.(|5l|) is automatically satisfied. Eq.(^) gives then the equation 

^xxxy ,xyyy — Oj 

whose solution is 

,xy ^ Fi{x + y) + F2{x - y). 

Therefore 

* = fi{x) + f2{v) + h{x -y) + U{x + y). 
The potential can then be written as 

V = vi{x) + V2{y) + V3,{x ~y) + V4,{x + y), 
where Vi = /f /4, i = 1, 2 and Vi = /2, i = 3, 4. Eq.(|5|) now gives 

«W4 + 2<ui) + 3v[v'4 - «U3 + wi) - Swiwg- 
(u2 W4 + 2<U2) - 3w2< + (v'^vs + 2w>2) - Sv'^v'^ = 0. 

This equation coincides with that reported in the review of Hietarinta p. Its most relevant solutions, 
together with the expressions of the invariant, are listed in table 2 below. 

In more general cases eq.(^) is not as easily solved. The standard approach is to restrict the attention 
to special cases where S{z) is either a homogeneous polynomial or a non-homogeneous one which is simply a 
power of the corresponding function in lower rank cases. Even with such restrictions, many solutions can be 
obtained. Moreover, in all the solutions found it has turned out that A{z) = Az^, with A a complex constant. 
In this way the coefficient 

zz + 2 5R{A(z)} 

of E in the Kahler potential (^9|) is always a Hermitian form in z. In the tables included in Section 6 are 
listed all the solutions found. 



5 Strong invariants generated by lower-order weak invariants 

In the present section we want to introduce an alternative technique to identify classes of strongly integrable 
systems. This approach works only in a restricted subclass of systems, but is useful to get insight into 
the structure of general solutions. The key point of the idea is based on the construction of a generic 
weak integrable system corresponding to the existence of a second-rank Killing tensor. Introducing the 
corresponding simplified Kahler potential in the general set of equations, we try to isolate single systems 
that satisfy them. The reader more interested in the general approach may wish to skip directly to section 

6 with the results. 



9 



5.1 Quadratic invariants 

The technique devised and apphed in paper |^ amounts, in short synthesis, in a conformal transformation of 
the form (p5|). In the second rank case, it turns out that it must be generated by an arbitrary holomorphic 
function S2{z) via the relation 

dz 



"'-''iTWy '''' 

In the new coordinates X, Y, the Jacobi potential G = E — V turns out to be of the form 

A{X) + B{Y) 
''^^^^^= \S,{X,Y)\ ' 

where A and B are arbitrary functions of their arguments, which can be thought of as the separation 
coordinates. The conformal factor appearing in ( |6^ can be written as 

|^2(^, Y)\ = ^Js2{w)S2{w) = F'{w)F'{w). (65) 

Moreover, we can write the transformation of the momenta in the form 

Px ^ Rpx + Q Py, ^gg^ 
PY = -Q Px + R Py, 

where 

Q = 

We can therefore write the expression of the second invariant in the Jacobi gauge, as 

Is{px,PY, X, Y) ^\{p\- pI) + B{Y) - A{X). (68) 

As remarked in the introduction, since to each Jacobi potential G pertains a specific dynamical system, we 
have that the found invariant is, in general, only a weak invariant for the standard Hamiltonian, in the sense 
that it provides a conserved quantity only at a given value of the energy (perhaps at zero energy). In fact, 
the Poisson bracket of the function ( |68| ) with the separated Hamiltonian 

Wx + pi) - A{X) - B{Y) , ^ 

Us = |„ ; (69) 

is 

{Is.ns} ^ ^,{PX \S2\,X~PY \S2\,y). (70) 

5.2 Quartic strong invariants generated by quadratic weak invariants 



In section 4, the system of equations (|5l|)"(|53|) was treated and many non trivial solutions were found. A 
general solution is however very difficult to find due to the complicated structure of the full system. Here 
we want instead to exploit the above results to generate new solutions simply starting from the assumption 
of the existence of a second-rank Killing tensor. 

As a foreword to this approach, let us make the following remark. Let us suppose to have made the choice 
of the S2 function generating the conformal transformation which gives the second invariant ( |68|) satisfying 
equation ([70|). Denoting by Tio the numerator appearing in the expression (|6^) of the Hamiltonian, the 
Poisson bracket of a generic higher-order invariant with T-Ls can be written as 



1 1~L 

\S2\ \S2? 

{IJ,n^}~E{IJ^S2\] ^ ' 



\S-2 
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In the case of a second-order invariant in the standard form of (|68|), this reduces to (|7^). On the other hand, 
let us introduce the null Hamiltonian 



n^f = {Hs - E)\S2\ =no - E\S2\. 
The Poisson bracket of the invariant with Ti._\f is 

{ij,n^} ^ {ij,no} - E{ij,\S2\}. 



(72) 



(73) 

Since 1521 is a non-nuU positive function everywhere, eqs.([7l|) and (^3|) are equivalent for what concernes 
the conservation of /,/. The conceptual difference relies on the fact that, whereas eq.(^3|) expresses the 
conservation of /,/ in the dynamics provided by Ti.^ at zero "energy", eq.([7l|) expresses the conservation of 
I J in the dynamics provided by Tis at arbitrary energy. Note that the true energy E enters into Ti.^/ as an 
arbitrary parameter. 

In the light of this argument, the procedure to look for new integrable systems is the following: regard 
the null Hamiltonian, written in the form 



1 



where 



GiX,Y) ^ A{X) + B{Y) + E\S2iX,Y)l 



(74) 



(75) 

as the system in which to find a new conserved quantity at zero energy. If the search of this new conserved 
quantity is successful, then it is a strongly conserved quantity for the original Hamiltonian Ti-s, which is 
therefore integrable at arbitrary energy. Going backwards with respect to the previous approach, we take 
the potential of eq.([75|) and use it to solve the system (^3[). Once we have the solution, we can write the 
integrability condition (|34|) in the form 



(76) 



Note that in the remaining part of this section we suppress the "tilde" over G and R to simplify the notation. 
The solution of this equation provides the explicit form of the functions A(w + w) and B(w — w) that allow 
for the integrability of the system. Finally, the integration of the equations for the trace (|3^) completes the 
solution. 

In the development of this approach, it can be useful to come back to the separating real variables. 
Equations (^) involving the Kahler potential therefore become: 



G(X,y) = i(/C^xx+/C,yy), 
RiX,Y) = i{lC^YY ~ IC,xx) - 
and the integrability condition ( [76|) can be written in the form 

mR})^XX - mR}).YY - 2{n{R})^XY 

3[((5{i?}),x - (5R{i?}),y)G,x - ((5R{i?}),x 



(77) 



i K, 



XY, 



G+ 



(3{i?}),y)G,y 



+ 



(78) 



Ki?}(G,xx-G,yy)-2 3?{i?}G 



XY 



= 0. 



The system 



finally becomes: 
4 



K,x = -^{mR}),x + mR}).Y))G- -(5R{i?}G,x + HR}G.y) 
K,Y - -\{mR}).x - (5R{i?}),y)G - ^(S{it:}G,x - 3fi{i?}G,y). 



(79) 
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To solve system (p3|), we need only consider the contribution to the potential (75) coming from the 
conformal factor: 

Gc E \S2{X,Y)\. (80) 

Recalling the relation ( |6^ ) between the conformal transformation and the S2 function in the quadratic case 
and the definition of the inverse transformation (54), the above source term can also be written as 



Gc=E S2{w)S2{w) = EF'{w)F\w). 
The solutions of system (^3|) for this term are: 

ICc = E {F{w)F{iu) + 2 ^{A{w)}) , 
Rc^~2E {F"{w)F{w) + A"(w)) 

so that the complete solution is therefore 



(81) 



(82) 



AdX 



dX 



BdY 



dY + E {F{w)F(w) + 2n{A{w)}) , 



i? = 2 {B{Y) - A{X)) ~2E {F"{w)F{w) + K" {w)) 



(83) 



This solution is subject to the integrability condition (|7^). In analogy with the general approach, substituting 
in it the solution ( p3| ) we get a second degree polynomial in E of the form (pO|). This equation must be 
satisfied for every value of the arbitrary parameter E^ so that each coefhcient must separately vanish. The 
three coefficients are in this case: 

A2 = -4i '^\^{F""F + A"")F'F' + 5F"'F"F'F + 3A"'F"F' + 2A"i^"'F'}, 

Ai = Ai 5{3(B,„ - A^^)F"F' - 3(^,„ + B,^){F"'F + A'") + 2{B - A)F"'F' 

- 2(A,„„ + B^^^){F"F + A") - (A + i?)F""F}, 

Aq = 2i 5{6(B „ - A,^){B^^ + A ,„) + 4(5 - A){A,^^ + B ^„)}. 

From the properties of the functions A and B, Aq is identically zero. Limiting ourselves to the simplest case 
A = 0, the condition that A2 vanishes turns out to be 



(84) 



+ 5F"'F")FF'} = 0, 

so that it is necessary and sufhcient that 

F""F' + 5F"'F" = 6aFF', a e I 

The general solution is 



(F 



f\4 



aF^ + + 6F - 



(85) 



(86) 



(87) 



where [3, 7, S, e are complex constants. We see that solution ( p7D is a subset of the general result (^8|) so 
that we can state that the simplified approach presented here consists in choosing a function in the set 

52 (z) = Vaz^ + iz^ + 5z + e. (88) 

The equation Ai = Q then picks the functions A and B that are compatible with strong integrability. The 
practical advantage with respect to the application of the general technique of section 4 is that equation 
^0 = is automatically satisfied. 
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5.3 Class iz"^ systems 

To illustrate the simplified approach, in this subsection we present a class of simple systems which are 
strongly integrable with a quartic invariant and, at zero energy, are separable with a quadratic invariant. 
Consider then the function 

S2{z) = (89) 
In the solution (^7|) it corresponds to the choice 



-1, 7 = ^ = 6 = 0. 



(90) 



This is the simplest polynomial form which does not satisfy the constraint (|60|) and therefore gives a po- 
tential which is not automatically integrable at arbitrary energy with a quadratic second invariant. The 
corresponding conformal transformation is given by 



Using polar coordinates, so that 



in terms of real variables, we get 



1 - i , 
w — H[z) — — — In z. 

v2 

z = X + iy = r cxp(i9) , 



r ~ \/ + y^, ^ = arctan(y/a;), 



X=-^(0 + lnr) 
r=-i=(^?-lnr) 



(91) 



(92) 



(93) 



We point out, in passing, that that in ||], in equations (102-104) corresponding to eq. (|93| ) above, there is a 
misprint and the coordinate r must be substituted with its natural logarithm. It follows that the potential 
given by 

^(6l + lnr)+B(6'-lnr) 



(94) 



is integrable at zero energy for arbitrary functions A and B. The common factor r ^ is due to the fact that, 
from the choice of eq. ( p9[) , we have 



|5(X,F)|=H-c 



_ „2 _ „\/2(x-y) 



The solutions are 



and 



/C = 4 



AdX 



dX + 4 



BdY 



dY + Ee^^''-^\ 



R = 2{{B{Y) - A{X)) -~2iE e^^^'^\ 
The integrability condition (^8|) becomes 

4 E e^^^-^) [A' - B" + 3V2{A' + B') + 4(yl - B)] = 0. 

This equation splits into 

A" + 3V2A' + 4A = 0, 
B" - 3V2B' + 4B = 0, 



with solutions: 



yl(X) = aie-^^+a2e-2^^, 
B(y) = 6ie^^ + 62e2^^. 



(95) 

(96) 
(97) 

(98) 

(99) 

(100) 
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The potential 

nr,g) = -"^^" + '^^'-°^^'":^^"" (101) 

ie therefore integrable. As a simple concrete example, let us take 

ai = 61 = 0, aa - &2 - -i, (102) 

so that the potential is 

26 I -20 

which is a known result mentioned in the review by Hietarinta. The second invariant has the form 

, , , , sh26l ch26l . 2-sh26' , , 

I=pi + ch29pl + 2—-prPe + i—rPe + 2 ■ (104) 

6 Classification of Hamiltonians admitting a quartic invariant at 
arbitrary energy 

In the tables included in this section we present all the solutions we have found of systems admitting an 
independent invariant quartic in the momenta. Some of them are new. For all cases we present the potential 
V (up to linear transformations of the coordinates and a linear rescaling of the potential itself) as well as the 
quartic invariant /, thereby making it possible to directly compare our results with Hietarinta's classification 
of 1987 For the cases which are superintegrable we also indicate, using the notation of j^, which of the 
real quadratic cases (1), (2), (4) or (7) a given system belongs to. 

A peculiar property of our solutions is that all but one have a conformal factor generated by function S of 
even degree in z, with real coefficients. As a consequence of this, the solutions obtained can be expressed in 
terms of the coordinate systems of separable potentials. In fact, since from (^) and ([s^), the transformation 
is 

H{z) = F-^{w{z)) = fdz [Siiz)]-^/^ (105) 



with S of the form 

S^iz) = ±{az^ + I3z + -f), aeM, /3,7eC, (106) 
we are actually taken back to the coordinate transformations to separable variables given by an S2 function 



satisfying the condition (60), for which the transformation is of the form 

H{z) = [ dz [S2{z)]-^/^. (107) 



As an exception to this situation, we have also found a solution with ^4 of third degree in z, by the method 
of the coupling constant metamorphosis (see Q) which is described in Section 6.2 below. 



6.1 Systems for which S{z) is of even degree 

The simplest case is that with Si of degree zero. This case corresponds simply to rotations and translations 
of the Cartesian coordinates. In table 1 we list the new solutions obtained and, for comparison, the quartic 
integrable potential of In table 2 the already known solutions of this class are listed together with 

the proper references. The first sequence of items correspond to the simple solution sketched at the end of 
section 4. All solutions are generated by S4 = ±1: actually, there is nothing fundamental in the choice of 
sign which has been made only for sake of comparison with already known results. 

We then have some solutions with 5*4 of degree two and four in z. In table 3, two new interesting classes 
of integrable systems are presented with ^4 of degree two. They are more conveniently treated using the 
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parabolic coordinates provided by the conformal transformation generated by iS'4 — 16z^. In table 5, a class 
of integrable systems generated by S4 — —z^ is presented, generalizing the separability in polar coordinates. 
Here, the choice of sign is dictated by the transformation used in the simplified approach of section 5.3 above, 
but to the same system one can also arrive with the general method using the "proper" transformation given 
by S'4 = z"*. 



6.2 A system for which S{z) is of degree three 

In order to describe how we get this new solution, we briefly recall the method based on the coupling constant 
metamorphosis We introduce a null Hamiltonian which, in analogy with eq.(|7^), is deflned as 

Um = [n - E)\S4\^ . (108) 

The second case of table 3 has a null Hamiltonian which, if expressed in the new coordinates, has the form 

1 



^Af - \{Vx +Py) - MX) - B{Y) + AC{X^ + Y'Y - AE{X' + Y^). 



(109) 



Then, if the term {X"^ + Y'^Y can be identified with the conformal factor associated to a new conformal 
transformation, we may interpret C as a new energy and as a new coupling constant. We get a new 
potential which is strongly integrable with an invariant of the same nature as before. This is the process 
known as coupling constant metamorphosis. 

In the present instance, we have that the transformation must be still of the form (105) or, equivalently. 



F'{w) = Sl'\ 

where, for simplicity, we still denote the new variable as w = X -\- iY, so that jwp 
conformal factor is assumed to be {X"^ + Y'^Y, we get 

\F'\^^\wf. 

Therefore, F{w) must be of the order w^, so that 

w^H{z) - z^^\ 

H'{z) ^ z-^l\ 
5*4(2) - 2^ 



X^ 



From ( |105| ) we get, 
so that, finally. 



(110) 
Since the 

(111) 

(112) 

(113) 
(114) 



We end up with a solution of third degree for S which had escaped a previous direct analysis. The new 
integrable Hamiltonian is 



2 J,/: ^ 



r2(l + cos|) r2(l-cos|) ' r3(l + cosf)3 ' r3(l-cos|) 



(115) 



where £* = C is the new energy and fj, = —E is the new coupling constant. Table 4 reports the form of the 
coordinate transformation and the expression of the second invariant. 
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S{z) = 1 

/C = 2Ex^ - |a;^ - 2bxy'^ - iy^ + 4c In \y\ 

V = Ax'^ + y'^ + bx + cy~'^ 

I=pt + 2pIpI + AVpl + 4byp^py + IQx'^pl + 2{2x'^ + y'^){h + Axf + 32ca;2y-2 
Superintegrable, case (4) and (7) 
S{z) = 1 

/C = 2Ex^ - 9y^/^{bx + c)-^ay^- ^dy^'^ - faa;^ 

V = ia(16a;2 + 9y^) + {bx + c)y-y^ + dy^^ 

I=pt + 2pIpI + AVpl + I2by^/^p^py + 32ax^pl + Max^V 
+l&bdx - 2560^0;* + ISft^yVs + lAAabxy'^/^ 

S{z) = -1 

/C = 2Ex^ - 2x^y^ - \axS - (f^ + 2/*) - if 2/' - + 4rfln|a;| 

V = x^ + &x^y^ + 8/ + a {^-y^ + x^y^ + ^ + 16^') + ^, 

I=pt+ Apl + 6x2y2 + ax^y + -^a^x^ + ^) ~ {^°-^^ + ^^^^^) 
+4a;^p2 ^ 4(2da;2 + a;^ + d^x''^ + ^ady + My^) + 

6 4 

= -1 

/C = 2^a;2 - 2x'^y'^ - ^y^ - i&cc^ - py^ + 4dln |x| - f j/"* + 4m In \y\ 

V = x^ + Qx^y^ + 8/ + bix^ + 4w^) + 4" + ^ + -ir> (Grammaticos et al. 1984) 

x/ x^ 

I = pI + 'ipl (^x^ + Qx^y"^ + ^ + ^) ~ '^^x^VP'^Pv + ^x^pI 

+4[l2cy2a;-4 + 2dx^ + x^ + 2bx^ + b'^x'^ + 2cx-'^ + 2cdx-^ + c^x'^^ + {d^ + 2bc)x-'^ 
+Ax^y'^ + 2mx^y-2 + 46x^y2 ^ 43,4^4 _^ 4^^2j 

Table 1: Systems for which S{z) is of zero degree. 
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7 Conclusions 



In this work we have explored the set of natural 2-dimensional Hamiltonian systems admitting a second 
invariant which is a polynomial in the momenta of degree four. The approach is based on the Jacobi 
geometrization described in where the linear and quadratic cases were also explored, and in [ p^ , where 
the cubic cases were treated. The approach allows a unified treatment of integrability at fixed and arbitrary 
energy, even if here we have limited ourselves to obtaining the list, as complete as possible, of strongly 
integrable systems. There are several areas in which this approach can be still useful. The most natural 
seems that of looking for systems admitting higher order invariants of which there exist only a few known 
examples. 

Moreover, other properties of invariants up to degree four are still worth of analysis. In fl^ we have 
examined the properties of weakly integrable systems with quadratic second invariants and found interesting 
behaviours both in the integrable regime and in the generic non-integrable one. One may wonder about the 
information higher order weak invariants can give on the dynamics of non-integrable systems. It is reasonable 
to believe that increasing the order can improve the reliability of these functions when interpreting them as 
effective invariants approximately conserved in regular portions of phase-space. It would also be interesting 
to look for other examples of quartic invariants which correspond to lower rank Killing tensors on fixed 
energy surfaces. 

References 

[1] V. N. Kolokol'tsov, Math. USSR Izvestiya 21, 291 (1983). 

[2] K. Rosquist and G. Pucacco, J. Phys. A 28, 3235 (1995). 

[3] S. Benenti, J. Math. Phys. 38, 6578 (1997). 

[4] V. I. Arnold, Mathematical methods of classical mechanics (Springer- Verlag, New York, USA, 1978). 

[5] C. Lanczos, The variational principles of mechanics (Dover, New York, USA, 1986). 

[6] G. Darboux, Archives Neerlandaises (ii) VI, 371 (1901). 

[7] L. S. Hall, Physica 8D, 90 (1983). 

[8] W. Sarlet, P. G. L. Leach and F. Cantrijn, Physica 17D, 87 (1985) 

[9] J. Hietarinta, Phys. Rep. 147, 87 (1987). 

[10] M. Karlovini and K. Rosquist, J. Math. Phys. 41, 370 (2000) 

[11] L. Samuelsson, Master's Thesis, Department of Physics, Stockholm University (1998), available at 
fittps/Zw-ww. fop. physto.se/exjobb/exjobb. htm, 

[12] V. V. Kozlov Symmetries. Topology and Resonances in Hamiltonian Mechanics (Berlin: Springer, 1991). 

[13] B. Grammaticos, B. Dorizzi, A. Ramani, J. Math. Phys. 25, 3470 (1984) 

[14] G. Pucacco and K. Rosquist, (2001), submitted 



17 



S{z) = 1 

/C = E{x^ + y2) - - f^{y) - f,{x -y)- U{x + y) 
V = vi{x) + V2{y) -\-vz{x - y) +V4:{x + y), 

ViiO = jfiiO (for i = 1,2), ViiO = yi'iO (for i = 3,4), (Hietarinta, 1987) 
I = pIpI + 2t^2P^ + 2{vi - V3)PxPy + ^VipI + h{x, y) 





= 0, U2 = ef , 


= e^ 


-y 


U4 = 0, /i = 2e^ + e2(^-2'). 




V\ 


= e~^, V2 = e", 


V3 = 




t;4 = 0, /i = 2e^ + e^f^-f) + 2e-f + 4e»-^. 




V\ 


= {), V2= e", V3 


= e^ 


5 


V4 = e-'^-y, h = 2e^ + e2(^-«) + 26"^ + e-2("^+^) - 2c-2y. 




V\ 


= 0, V2= e", V3 


= e^ 


-y 

5 


t;4 = e-(^+2')/2, /i = 2e^ + e^^^-'') + e-(^+^') - 2e(^-32')/2. 




V\ 


= e~^^, V2 = e^, 


t^s = 




-y, V4 = 0, h = 2e^ + e^^^-'') + Ae-'^^'+y. 




V\ 


a a 


U3 = 






16a6 


(x 


— y)^' (x + j/)^ ' x^j/^ (x^ — y^)** (x^ — y^)'^ 


(x^ — y^)^ 



+8^77^ + (32am + 72c^)x\^/^ + 8nmx"\^/^ + 72anx"\^ 

+4;r.r- 1 + :-!2«r.r-(/--''''(9(/- + .r-) + :V2(r.r-(':}tr + 'i.r^) 



2 ' ^ ' V 3 

I =Pt+^Pli 2°'^^ ~^ '^^^y + + ^ ) ~ -^dx^PxPv 



2dh2yx~'^ + 2hih2X-~^ + /lix"^^ + {hj + a/i2)x"^ + ia^x'* - :^c?^a:'^ 

4 18 



+4 

-^(fx^y"^ - ^adx'^y + ^dhiy 



S{z) = 1 

K. = 2Ex'^ - 9y4/3(cx2 + d+ ^y^^^) - a^x^ + 3y^) + 4n In |x| 

9 Tl 

V = -cy"^^^ + (cx^ + d)y~^^^ + my'^/^ + a{9y^ + Ax"^) + — , (Grammaticos et al. 1984) 
2 x^ 

I = p'i + 2pIpI + 4p| [{cx^ + d)y-^/^ + my"^^^ + a{9y^ + \x^) + ^) + 2Acxy^l'^p^py 

+4plUax'^ + nx~^) + IQcmx^ + 32adx'^y^^/^ + 8 
cn 



S{z) = -1 

/C = 2£;x2 - a (ix^ + §/) - d (Ix^y + §y^) + Aht In |x| - ^ 
V = \a{x^ + 162/2) + d i^y^ + x^y] + % + ^, (Grammaticos et al. 1984) 



Table 2: Systems for which S{z) is of zero degree (continued). 
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S{z) = 16z^ 
/C = |£(2a;2 + y2) + /(X) + g{Y) + C{XYY 



/(X) = -36 2-l/3cX4/3 _ 182l^aX8/3, = -36 2-1/3^^4/3 _ 18|l^^y8/3^ „ ^ 4/3, 

+ c(r + X 

)-i/3+rf(^_a,)-i/3 ^ _ 21/3(7 



,2/3 ' 



A; 



36 



7 = J2 - 8A;t/i/3 [p^pe - ^{a{r + xfl"^ + 6(r - xf/"^ + c(r + x)-^l^ + d{r - x)-^^)] 
-8k [a{r - x)2/3 + b{r + x)2/3] , 



J = PyPe + 
Pe = xpy - ypx- 



{r+x)[b{r-x)^^'^+dir-x)-'-^'^]-{r-x)[a{r+xy/'^+c{r+x)-'-^'^] 



+ 2 



kx 



f{X) = ^CX^ + 8alogX - ^cX-\ g{Y) ^ ^CY^ + SblogY - ^dY-\ a = 4 



r 



1. 



C 



-kr , k = — — , 



r + X r — x {r + .t)3 (r — .t)3 
1 = 3"^- ky^ Ip^ps - ^ f ^ + + + r-^ 

I'xfu ^ \r+x r—x {r-\-xy (r— x) 

-k [o(r - xf + 6(r + xf] , 
J = PyPe + i [(r + a;)[^ + - (r - a:)[^ + 



+ ^kxy'^. 



Table 3: Systems for which S{z) is of second degree. 
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Siz) = 64z^ 

K = ^E{Ax^ + 3?/) + f{X) + g(Y) - ^^l{X^ + Y^) 

X = V2rV4cos| = ri/yi + cosf, Y = V^r^/^sinl = r^'^^l-cos^, 

r = + y2 = (X2 + y2)2/4, 6» = arctan(y/a;). 

f{X) = 32a \ogX-lcX-\ 

giY) = 32b logF-fdF-^. 

H a b c d 

" r ^ r2(l + cosf) ^ r2(l-cosf) ^ r3(l + cos|)3 ^ r3(l-cosf)3' 

7 = j2 + rin-r2, 

J = A + B(y) - A{x) - n{x^ - y2)3 + 4^(x2 _ y2^^ 

^ = 7? [^°^ f ('■^^'r - Pe) - sin ^rprPe] , 
^=7? [sin |(r2p2 _ p2) ^ ^rprPe] , 
Ti = 16 W x^ys^ 

r2 = 16 w [n {2X^Y^ + sx^v^ix^ + ^^)) - s/ix^y^ - ay^ - ftx^]. 

H is the Hamiltonian of the system, i.e. H = ^{p^ + ^) + V, 
P0 = xpy - yp^, 
rpr = xpx + ypy. 

Table 4; Systems for which S(z) is of third degree. 



S{z) = -z^ 

JC = Er^ + f{X)+g{Y) 

X=^{e + \nr), Y=^{0-lnr). 

A{X) = aie-^^ + a2e-2^^, A = f" /A 

B{Y) = bie^^ + 626^^^, B = g"/A. 
y _ae^ + be~^ _^ ce2^ + de~2^ 

I=pt + (ce2« + de-2«)p2 + 2 [ae^ - be"^ + i(ce2^ - de-^^)] PrPe 
+ (ae^ + be-') + ^{ce^' + de-^')] pj 

Table 5: Systems for which S{z) is of fourth degree. 
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